AN EQUIVALENCE RELATION ON WAVELETS IN HIGHER 
DIMENSIONS ASSOCIATED WITH MATRIX DILATIONS 
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Abstract. We introduce an equivalence relation on the set of single wavelets of 
L 2 (M™) associated with an arbitrary dilation matrix. The corresponding equivalence 
classes are characterized in terms of the support of the Fourier transform of wavelets 
and it is shown that each of these classes is non-empty. 



1. Introduction 

Let n > 1 be an integer. An n x n matrix A will be called a dilation matrix for R n if 
A preserves the lattice Z n , that is, ATP C T n , and all eigenvalues of A have absolute 
value greater than 1. These conditions imply that A has integer entries and | det A\ is 
an integer greater than 1. Let B = A*, the transpose of A, and a = | det A\ — | det B\. 

Since ATI 1 is a normal subgroup of the abelian group T n , we can form the cosets of 
ATP in T n . It is a well-known fact that the number of distinct cosets of ATP in TP is 
equal to a = | det A\ ([6], [11])- A subset of TP which consists of exactly one element 
from each of these a cosets will be called a set of digits for A. It is easy to see that if 
K is a set of digits for A, then so is K — fi, where fi G K. Therefore, we can assume, 
without loss of generality, that G K. 

A finite set of functions \1/ = {ip 1 , ip 2 , . . . , ip L } C L 2 (M. n ) is said to be a multiwavelet 
of L 2 (M. n ) associated with the dilation A if the system 

{^. k ■= a? l2 ^\A 3 ■ -k) : j e Z, k G Z n , 1 < I < L} 

forms an orthonormal basis for L 2 (M n ). If there is a single function ip G L 2 (R n ) such 
that 

{iP jk ■= a j/2 ijj(A j ■ -k) : j eT,k G T n } 

is an orthonormal basis for L 2 (MP), then we say ip to be an A-wavelet. 

A. Calogero [4] provided a characterization of multiwavelets associated with a gen- 
eral dilation matrix. This result extends to the higher dimensional situation a well- 
known theorem for dyadic wavelets in L 2 (M), which was proved independently by 
G. Gripenberg and X. Wang (see Chapter 7 of [7] and reference cited there). Dai, 
Larson, and Speegle ([5]) proved that, associated with every dilation matrix, there 
exist wavelets ip of L 2 (MP) of the form = xk, where K is a measurable subset of 
MP. Such wavelets are called MSF (minimally supported frequency) A- wavelets and 
we shall call the corresponding set K an A-wavelet set. The theorem of Calogero for 
single wavelets of L 2 (MP) is the following. 
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Theorem 1.1. A function ip G L 2 (IR n ) is an A-wavelet for L 2 (IR n ) if and only if the 
following are satisfied: 

(i) IHh = l, 

(ii) EI^'0I 2 = 1 fora.e.£eW\ 
jez 

(iii) HB j OH Bj (t + 2gvr)) = for a.e. £ G R n , q E Z n \ BZ n . 

j>o 

We use the following definition of Fourier transform: 

f(0 = [ f(x)e-^dx, £ G R n . 

The A-wavelet sets (equivalently, the MSF A-wavelets) admit a simple characteri- 
zation in terms of two geometric conditions. 

Theorem 1.2. A set K C R n is an A-wavelet set if and only if the following two 
conditions hold: 

(i) {K + 2kix : k G Z™} is a partition ofR n , and 

(ii) {B^K : j G Z} is a partition ofM. n . 

Given an A-wavelet ip, we define a sequence of closed subspaces of L 2 (M. n ) as follows: 
Vj = span{V/, fc :l<j,ke Z n }, j G Z. 
It is easy to verify that these subspaces satisfy the following properties: 

(i) Vj C V j+1 for all j G Z, 

(ii) / G Vj if and only if f(A-) G V j+1 for all j G Z, 

(iii) Uj- ez ^ is dense in L 2 (R n ), 

(iv) rij e z^ = {0}, and 

(v) Vo is invariant under translation by elements of the group Z n . 

In view of the last property, it is natural to ask the question of the existence of 
other (larger) groups of translations under which the space Vo remains invariant. For 
each non-negative integer r, we shall prove the existence of A-wavelets for which the 
corresponding space Vo is invariant under translation by elements of the form A~ r k 
for all k G Z n . 

For y G R n , let T y be the (unitary) translation operator on L 2 (R n ) defined by 
Tyf{x) = f(x — y). We consider the following groups of translation operators: 

Qr = {T A -r k ■ k G Z n }, r G N U {0}, and = {T y : y e R n }. 

Let Q be a set of bounded linear operators on L 2 (R n ) and V a closed subspace of 
L 2 (R n ). We say that V is ^-invariant if Tf G V for every / G V and T G Q. 

Let L r denote the collection of all A-wavelets such that the corresponding space V 
is ^-invariant. It is clear that L is the set of all A-wavelets, and 

Lo Z) Li Z) L2 Z) • • • D L r Z> L r _|_i Z) • • • Z) L^. 

These inclusions naturally give rise to an equivalence relation on the collection 
of A-wavelets. The equivalence classes are given by Ai r = L r \ L r +i, f > 0, and 
A^oo = Lqo. Therefore, A4 r , r > 0, consists of those A-wavelets for which Vo is 
^-invariant but not ^ r+ i-invariant. 
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This equivalence relation was defined in [10] for dyadic wavelets in L 2 (M), where the 
author characterized the equivalence classes. He also proved that M r , r — 0, 1,2,3, 
are non-empty. Later, in [2], [9], examples of wavelets of L 2 (IR) were constructed for 
each of these equivalence classes, by different methods. This equivalence relation was 
extended by the author to the wavelets of the Hardy space H 2 {R) in [1], where it was 
also shown that the corresponding equivalence classes are non-empty. 

In section 2 we characterize the equivalence classes in terms of the support of the 
Fourier transform of wavelets and in section 3 we prove that each of them is non- 
empty. 

2. Characterization of the equivalence classes 

For j G Z and y G M™, define the dilation operators Dj and translation operators 
T y on L 2 (R n ) as follows: 

D j f(x)=a?l 2 f(Aix) and T y f(x) = f(x-y). 

In this notation we have ipj k = DjT k vp. Observe that for j G Z and y G W l , we have 

TyDjfix) = Djfix-y) 

= a^ 2 f{A^x-y)) 

= '/;, ,//' 2 /i.V'.r) 

= P/IajJU) 

Hence, we have the following commutation relation: 

(2.1) T y Dj = DjT A j y , j' G Z and y G R n . 

We define another sequence, {Wj : j G Z}, of closed subspaces of L 2 (IR n ) by 
Wj = Vj + i Q Vj, the orthogonal complement of Vj in V j+1 . Clearly, 

(2.2) I 2 (l") = K„®(0f 3 ). 

i>o 

Observe that W = span{'0(- — k) : k G Z n }. Therefore, the wavelet ip belongs to the 
space Wo- For this reason, first we characterize this space. 

If / G Wq, then there exists a sequence {c k : k G Z n } G / 2 (Z n ) such that / = 
Ylk& n c fcV ; (' — k). Taking Fourier transform of both sides, we get /(£) = ?n(£)V>(£) 
for some 27rZ n -periodic function m G L 2 (T n ), where T n = [— n,7r] n . Conversely, if / 
is of this form, then by considering the Fourier series of m and reversing the above 
steps, we see that / G Wq. Therefore, we have the following characterization of Wq. 

(2.3) W = {/ G L 2 (W l ) : /(£) = m(f)^(0 for some m G L 2 (Y n )}. 
We begin with a few lemmas. 

Lemma 2.1. Fix a non-negative integer r , and let p, q be integers with p < q. Then 
the space V p (resp. W p ) is Q r -invariant if and only ifV g (resp. W q ) is Q r+q _ p -invariant. 

Proof. By definition, / G V p if and only if D q _ p G V q . Suppose that g G V q and find 
f E V p such that D q - P f = g. Now, from the commutation relation (2.1), we have for 
k G Z n 

T A -(r+ q -v) k D q -pf(x) = D q - p T A -r k f(x). 
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The lemma now follows from this equation. □ 
The structure of the space V can be quite complicated; in some cases it may not 
even be generated by the Z n -translates of a finite number of functions. On the other 
hand, Wq is generated by the Z n -translates of a single function, namely, the wavelet 
ip. The following lemma allows us to work with Wq instead of Vq. 

Lemma 2.2. The space Vq is Q r -invariant if and only ifWo is Q r -invariant. 

Proof. Suppose that Vq is ^-invariant. By Lemma 2.1, V\ is (? r+ i-invariant. In 
particular, V\ is ^-invariant. But V± — V © Wq. Hence, Wq is ^-invariant. 

Conversely, let Wq is ^-invariant. Then, again by Lemma 2.1, Wj is Q r ^-invariant 
for all j > 0. In particular, Wj is ^-invariant, hence so is ®j>oWj. The £ r -invariance 
of Vq now follows from the decomposition (2.2). □ 

Lemma 2.3. Let M be a dilation matrix and Km be a set of digits for M . Then 

e" 27 ™^ 1 ^ = | det M\5 0l/ , for v e K M t, 

where Km* is a set of digits for the dilation M l . 

Proof. This is the orthogonal relation for the characters of the finite group Z"/MZ" 
(see [8]). Observe that the mapping 

/i + MZ n h-> e -^W~ W> ) v e KMt 

is a character of the (finite) coset group Z n /MZ n . If v — 0, then there is nothing to 
prove. Suppose that v ^ 0, then there exists a // G K M such that e ~ 2m ( M M ' ^ ^ 1. 
Since K M is a set of digits for M, so is K M — //. Hence, 

(2.4) e- 2 ™^ 1 ' 1 ^-^'^ = 

Now 

= e -27rj(M-V'^> . ^ e "27rj(M-V,^> ) (2.4). 

Therefore, 

J2 e -2-<M-W> = o, since e - 2 ™< M ^>> + 1. 
/j.eK M 

□ 

Applying Lemma 2.3 to the dilation matrix A r , we get the following simple but 
useful fact. 

Lemma 2.4. Fix r e N and let I e Z n \ B r Z n . Then there exists k e Z n swc/i t/iat 
(B- r l,k)^Z. 

Proof. Since A is a dilation matrix, so is A r . From Lemma 2.3, we have 

e^™^^ = 0, veKsr\ {0}. 



e -27rj(M-V^>_ 
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If I G Z n \ B r Z n , then 1 = B r p + u' , where p e Z n and v' G ^ \ {0}. Therefore, for 
pi G K A r, we have (A~ r /i, I) = (A~ r fi, B r p + v') = (fi,p) + (A~ r fi, i/); hence, 

(2.5) ^ e~ 2ni{A ~ r = ^ e~ 27ri<A ~' > ' l/> = 0. 

Now, if (B~ r l, k) = (A- r k, I) e Z for all k G Z n , then each term in the sum in (2.5) 
would be equal to 1, a contradiction. □ 

Let U and V be two measurable subsets of M™. Then U is said to be 27rZ n - 
translation equivalent to V if there exists a measurable partition {Uk '■ k G Z n } of U 
such that {C/fe + 2kn : /c G Z n } forms a partition of V . A function / on W 1 is said to 
be 2nZ n -periodic if /(x + 2mvr) = /(x) for a.e. x G M n and all m G Z n . 

Let ^ be an A-wavelet. Denote E = supp t[j. For fc G Z n , define 

£(Jfe) = {£ G £ : f + 2kn G £} = £ n (E + 2A;vr). 

Theorem 2.5. Lei r G N and ip be an A-wavelet. Then ip G L r if and only if 
E(k) = for every keZ n \ B r Z n . 

Proof. Necessity. Let ip G L r . That is, Vq is ^-invariant, hence by Lemma 2.2, Wq 
is ^-invariant. Since ip G Wq, the £ r -invariance of Wq implies that — A~ r k) G Wq 
for all k G Z n . Taking Fourier transform, we get e~*' r '' fc ^ G Wo for all fc G Z n . 
By the characterization of Wq (see (2.3)), this is equivalent to saying that for each 
k G Z n , there exists a 27rZ n -periodic function /i fc G L 2 (T n ) such that e~^ B r ^' fc V(£) = 
^(0^(0- Note tl:L at if £ G -B (so that V>(£) 7^ 0), then we have 

(2.6) e- i < B_r ^>=^(0- 

Now suppose that there exists / G Z n \ 5 r Z n such that i?(Z) 7^ 0. Hence, there exists 
a non-trivial set F such that both F and F + 2/7r are subsets of E. Now, for a.e. 
(Gf and for all k G Z n , we have 

e -i<B-K+2br),fc>^ + 2/7r) = ^ + 2/7r) ^ (e + 2/7r) 
=}► e -i(B-^k) e -2ni(B-ri, k )^ + 2/?r) = ^(f + 2/vr) 

(since is 27rZ n — periodic) 
=> e -z«HB- r i,k) = 1 ( by (2.6), and since £ + 2ln G f?) 
=^ (B~ r l,k)eZ forallA;GZ n . 

This is a contradiction to Lemma 2.4. 

Sufficiency. Suppose that E(k) = for all k G Z n \ 5 r Z n . To establish ^ G L r , 
it is sufficient to show that for all k G Z n , we have 

(2.7) e -*< B -€.*>^(0 = ^(0^(0 

for some 27rZ"-periodic function v k G L 2 (T n ). Define the map 

r : R n — > T n , r(x) = x + 2m(x)7r, 

where m(x) is the unique element in Z n such that x + 2m(x)7r G T™. Since the relation 
2~2mez n IVKC + 2m7r)| 2 = 1 a.e. holds for every A-wavelet, it is clear that r : i£ — > T™ 
is onto. First, we want to find a subset H of E such that r : H — > T n is a bijection, 
that is, H is 27rZ n -translation equivalent to T n . 
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Fix f G T n . Define Z(g) = {k G Z n : £ + 2kn G E}. Notice that for a.e. 
£, Z(£) is non-empty. For £ G E, we choose fc^ = 0. Now, let £ ^ E. If there 
is an I = l 2 , . . . ,l n ) G Z(£) such that li > 0, then we choose = k, where 
ki = infj/x : / = l 2 , . . . , l n ) G Z(£) and l\ > 0}. Otherwise we choose k% = k, 
where k\ = supj/x : I = (li, l 2 , . . . , l n ) G Z(£) and /i < 0}. Now define the set 

(2.8) H = {S + 2k ( p:Sel n }. 

Observe that H is 27rZ ra -translation equivalent to T n , by construction. For k G Z n , 
we define the function on H by 

(2.9) u k (0 = e-^"^) fcefl) 

and extend 27rZ n -periodically to the whole of M n so that (2.7) holds clearly on H. To 
complete the proof we need to show that (2.7) also holds on E \ H. 

Now, for almost every £ G E \ H, there exists (' 6 and ^ G Z n such that 
£ = f + 2/ ? 7T, since if is 27rZ™-translation equivalent to T n . That is, E(l ( ) ^ 0. 
Hence, by hypothesis, 1% = B r p^ for some p$ G Z™. Therefore, we have 

= e-^" r «' fc ^(£' + 2^7r) 

= + 2Z € tt), (by (2.9), since £' G if) 

= ^(£-2/^(0 

— ^(0^(0 (since z/ fc is 27rZ n — periodic). 

This completes the proof of the theorem. □ 
Theorem 2.5 allows us to characterize the equivalence classes in terms of the support 
of the Fourier transform of the wavelets. 

Theorem 2.6. (a) The equivalence class M.00 is precisely the collection of all 
MSF A-wavelets. 

(b) An A-wavelet ip G M r , r > I, if and only if E(k) = for all k G Z" \ B r Z n 
but there exists I G Z" \ B r+1 Z n such that E{1) ^ 0. 

(c) An A-wavelet tp G A4 if and only if there exists I G Z n \ BZ n such that 

Proof. Item (a) is an easy generalization of the corresponding result for the one- 
dimensional dyadic wavelets, which is proved in [10]. Item (b) follows from Theo- 
rem 2.5. It also follows from Theorem 2.5 that an A-wavelet ip G La if and only if 
E(l) = for every I G Z n \ BZ n . Since L consists of all A-wavelets, (c) follows from 
this and the fact that M — L \L 1 . □ 

3. Construction of wavelets in M r 

Let D be a subset of M™ satisfying (ii) in Theorem 1.2. For E C M n , we define the 
translation projection r and dilation projection d as follows: 

t(E)= [J {(E + 2kn) nT n },and 
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d(E) = \J{(BiE)nD}. 

The following result gives sufficient conditions for a set S C W l to be a subset of 
an A- wavelet set. We refer to [3] for the proof. 

Theorem 3.1. Let A be a dilation matrix, B = A 1 , and S C M. n . Then S is a subset 
of an A-wavelet set if the following conditions hold: 

(a) S n (S + 2k-n) = for all k G Z n , 

(b) (B'S) n S = for all j G Z, 

(c) there exists e > such that N e (0) C\t(S) = 0, where N € (x) is the ball of radius 
e around x, 

(d) D \ d(S) has a non-empty interior. 

We shall also require another result proved in [3] for expansive matrices, i.e., n x n 
matrices such that all eigenvalues have modulus greater than 1. Hence, a dilation 
matrix is an expansive matrix which preserves the lattice Z n . 

Lemma 3.2. Suppose A is a dilation matrix, k G Z n \ {0} 7 and p G Z \ {0}. For 

a G Z n and j G Z ; let f a (x) = B~ p x — 2an and gj(x) = B~ix — 2kn. Let Y be the 
collection of all fixed points of T U Q , where 

(3.1) T = {f a : a G Z"} and Q = { Qj : j G Z}. 
Then for any y (jLY , there exists e > such that 

(3.2) fa(N e (y)) D N e (y) = 9j (N e (y)) n N e (y) = for all a G Z n ,j G Z. 

Using Lemma 3.2, we show the existence of a set I C M™ satisfying some properties, 
which will be crucial for the construction of wavelets in M. r . 

Proposition 3.3. Let k G Z n \ {0} and p G Z \ {0}. TTiere exists a measurable set 
I dW 1 satisfying the following properties: 



(3.3) |/| > 0, 

(3.4) t(I) n t(B- p I) = 0, 

(3.5) d(I)f]d(I + 2B p kn) = 

(3.6) / U (B~ P I + 2kn) is contained in an A-wavelet set. 



Proof. Let Y be the collection of fixed points of T U Q as defined in Lemma 3.2. 
Notice that the fixed points of T is precisely the lattice 2n(B~ p — J)^ 1 Z ra and those 
of Q is the set {2tt(B~^ - I^k : j G Z \ {0}}, which has limit points {0, -2kn}. 
Choose any y £ Yu (2vrZ n ) U (2nB p Z n ). By Lemma 3.2, there exists an e > such 
that (3.2) holds. That is, if we take / = N e (y), then I satisfies (3.3), (3.4), and 
(3.5). To prove (3.6) we shall show that S — IU (B~ P I + 2/c7r) satisfies the conditions 
(a)-(d) of Theorem 3.1. Equations (3.4) and (3.5) clearly imply (a) and (b). Since 
y G" (27rZ n ) U (2nB p Z n ) ) by chosing e' < e sufficiently small, if necessary, and taking 
/ = N' e (y), we observe that (c) and (d) are also satisfied. Hence, there is an A-wavelet 
set containing / U (B~ P I + 2kn). □ 
Now fix an integer r > 0. Choose k — k r G B r Z n \ B r+1 Z n and p — 1. Then by 
Proposition 3.3, there is a set I C M. n of positive measure and an A-wavelet set W 
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such that 

(3.7) r(/)nr(fi- 1 /) =0, 

(3.8) d(I) n d(I + 2Bk r ix) = 0, 

(3.9) / U (B^I + 2k r 7i) C W. 

Let J = W \ (I U (f? _1 f + 2k r 7i)), and define the function ?/v as follows: 

' ^ if £ e f U (5- 1 /) U (fi- 1 / + 2A; r vr) 
if £ G f + 2Bk r TC 
1 if £ G if 
otherwise. 



(3.io) Vv(0 



Theorem 3.4. For eac/j integer r > 0, an A-wavelet belonging to the equivalence 
class M. r . 

Proof. To prove that ip r is an A-wavelet, we verify that ip r satisfies (i)-(iii) of Theo- 
rem 1.1. Condition (i) is verified once we show that X^mez™ lVv(£ + 2mir)\ 2 = 1 a.e. 
For, 



r \\2 



[ \MO\ 2 d£= I I^(e + 2m7r)| 2 ^=(2vrr, 



hence, ||Vv||i = 1- Let p(£) = XlmeZ" IV't-C^ + 2kn)\ 2 . Since p is 27rZ n -periodic, it 
suffices to prove that p(£) = 1 a.e. on the wavelet set W which is 27rZ n -translation 
equivalent to T\ Notice that 

W = J U I U (B^I + 2A; r 7r) and supp Vv = W U (B -1 / ) U (f + 2Bkrir). 

If £ G J, then £ + 2mi{ G supp ?/v if and only if m — 0, hence p(£) = 1. If £ G f , 
then £ + 2m7r G supp if and only if m = 0, or m = Bk r . Hence, p(£) = |Vv(£) 1 2 + 
|Vv(£ + 2Bk r n)\ 2 = (^) 2 + (-^) 2 = 1. Finally, if £ G B~ l I + 2k r n, then £ + 2mn G 
supp ^ if and only if m = 0, or m = —k r so that p(£) = |?/v(OI 2 + l^(C — 2/c r 7r)| 2 = 

<7!) 2 + (A) 2 = 1 - 

The proof of (ii) is similar. 

Now we prove (iii). The term ^j r (B^)^j r (B^ (£ + 2qn)) is non-zero only when both 
B^ and + 2qn) are in the support of ip r . Since -B J (£ + 2qn) = fP£ + 2B^qii, it 
is clear from (3.7), (3.9), and the definition of ip r that either B^q = k r or B*q = Bk r . 
In the first case, we get that f? J '(£ + 2qir) G I + 2Bk r ir and fP£ G f if and only if 
B j ~ l (l; + 2gvr) G £ _1 J + 2A; r 7r and fi^ 1 ^ G B~ X I . Hence, the sum in Theorem 1.1 
(iii) is equal to (^=) (—75) + ("75) ("75) = 0- The second case is similar. 

Therefore, by Theorem 1.1, it follows that ip r is an A-wavelet. 

Finally, we have to show that ipr G M. r . Recall that E(k) = supp ip r D (supp ip r + 
2kn). From the definition of ip r , it is clear that E(k) is non-empty if and only if 
k = 0,±k r ,±Bk r . By our choice, k r G B r Z n \ B r+1 Z n . Hence, by Theorem 2.6(b) 
and (c), ip r G M r , r > 0. □ 

As we mentioned earlier, it was proved in [5] that there are MSF A-wavelets of 
L 2 (IR n ) associated with every dilation matrix. Hence, A4oo is non-empty. We just 



AN EQUIVALENCE RELATION ON WAVELETS IN HIGHER DIMENSIONS 



9 



proved that Ai r is non-empty for r > 0. Therefore, each of the equivalence classes of 
A- wavelets is non-empty. 
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